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^ [ Abstract 
ON , 

■ The di-photon decay channel of the hghtest Higgs boson is considerd as a probe to 

j. , explore CP violation in the Minimal Supersymmetric Standard Model (MSSM). The 

I scalar /pseudo-scalar mixing is considered along with CP violation entering through 

' the Higgs-sfermion-sfermion couplings, with and without light sparticles. The impact 

. of a light stop on the decay width and Branching Ratio (BR) is established through 

^P. I a detailed study of the amplitude of the process Hi ^77. The other sparticles have 

' little influence even when they are light. With a suitable combination of other MSSM 

^ ■ parameters, a light stop can change the BR by more than 50% with a CP-violating 

^ , phase (j)^ ~ 90°, while the change is almost nil with a heavy stop. 

1 Introduction 

Despite its success to describe the physics of elementary particles there are strong hints that 
the Standard Model (SM) is only an effective theory valid up to the TeV range, and new 
physics is needed to explain particle dynamics (much) beyond this energy scale. Among 
many others Supersymmetry (SUSY) is one of the most favoured scenarios and will be 
explored for in all possible ways at the upcoming Large Hadron Collider (LHC) at CERN. 
Furthermore, one of the main tasks of the LHC is the search for Higgs bosons, i.e. the 
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determination of the underlying mechanism of Electro- Weak Symmetry Breaking (EWSB). 
Precision measurements at CERN LEP and SLAC SLC prefer a light Higgs particle, which 
is indeed predicted in SUSY models. The scalar potential of the MSSM conserves CP at tree 
level [1], because SUSY imposes an additional (holomorphic) symmetry on the Higgs sector 
of a general two-Higgs doublet model, that enforces flavour conservation in tree-level neutral 
currents and absence of CP- violating scalar/pseudo-scalar mixing in the Born approximation. 
Beyond tree-level, the CP invariance of the Higgs potential may in principle be spontaneously 
broken by radiative corrections when the Vacuum Expectation Values (VEVs) of the two 
Higgs doublets develop a relative phase [2, 3]. According to the Georgi-Pais theorem [4] 
though, this type of CP violation requires a very light Higgs state, which is essentially ruled 
out by experiment [5]. 

On the other hand many new MSSM parameters can well be complex and thus explicitly 
break CP invariance. Beyond the Born approximation these new CP-violating phases induce 
CP violation also in the Higgs sector [3, 6]. The possibly complex parameters include: (i) 
the higgsino mass term /i, (ii) the soft SUSY-breaking gaugino masses Ma (a = 1,2,3), 
(iii) the soft bi-linear term Bji and (iv) the soft trilinear Yukawa couplings ^/ of the Higgs 
particles to scalar fermions of flavour /. In general, all of these new phases are independent. 
However, when imposing universality conditions at a unification scale Mqut all the three 
gaugino masses have a common phase as well as all the trilinear couplings Af have 
another common phase, i.e., four independent phases remain: those of B/x, Ma and Af. 
Furthermore, the two U{1) symmetries of the conformal-invariant part of the MSSM may 
be employed to re-phase one of the Higgs doublet fields and the gaugino fields such that 
Ma and are real [6, 7]. In this paper we will work within this setup with two independent 
physical phases, which we take to be arg(/x) = 0^ and axg{Af) — (j)Af 

The new CP-violating phases in the MSSM are severely constrained by bounds on the 
Electric Dipole Moments (EDMs) of electron, neutron and the Hg atom. In order to avoid 
problems with phases associated with the sfermions of the first and second generations one 
may deviate from exact universality and consider to be diagonal in flavour space with 
vanishing flrst and second generation couplings [8]. In general, the constraints are rather 
model dependent and there have been several suggestions [9]-[14] to evade these constraints 
allowing large CP-violating phases of 0{1). One possibility is to arrange for partial cancel- 
lations among various contributions to the EDMs [10]. In this scenario, it has recently been 
pointed out that for large trilinear scalar couplings Af, phases (p^ ~ 0{1) can be compatible 
with the EDM bounds [11]. Another option is to make the flrst two generations of scalar 
fermions rather heavy, of order a few TeV, so that the one-loop EDM constraints are auto- 
matically evaded, however, two-loop contributions of third generation scalar fermions may 
still be large [8] . (A detailed analysis of the so-called CPX scenario with very heavy flrst and 
second generations squarks is available in Ref. [12].) As a matter of fact, one can consider 
so-called effective SUSY models [13] where decoupling of the flrst and second generation 
sfermions are invoked to solve the SUSY Flavour Changing Neutral Current (FCNC) and 
CP problems without spoiling the naturalness condition. Furthermore, the restrictions on 
the phases may also disappear if lepton flavour violating terms in the MSSM Lagrangian are 
included [14]. In conclusion, this means that large phases cannot be ruled out and therefore 



2 



we analyse the full range 0° < (0^, 0a) < 180° in the following. 

The CP- violating SUSY phases can in principle be determined directly in SUSY particles 
production and decay at high energy colliders [6], [15] -[23] or indirectly via their radiative 
effects on the Higgs sector [6, 17]. Here we focus on the di-photon decay mode. Hi 77, of 
the lightest neutral Higgs boson Hi, which involves direct, i.e leading, effects of the SUSY 
phases through couplings of the Hi to SUSY particles in the loops (sec Fig. 1) as well as 
indirect, i.e. sub-leading, effects through the scalar/pseudo-scalar mixing yielding the Higgs 
mass-eigenstate Hi. The di-photon decay mode is important for the study of CP- violating 
effects in the MSSM Higgs sector for two reasons. On the one hand, it is the most promising 
channel for the discovery of a light neutral Higgs state with mass between 80 and 130 GeV 
at the LHC [24, 25]. On the other hand, the coupling strength of the dominant CP- violating 
terms of the di-photon decay width which depend on /x and Af (with / = b,t,T, hereafter) 
is of the same order, 0{a^), as that of the CP-conserving ones. 

The entire gg/qq //i ^ 77 process can be factorised exactly into three parts: the 
production process, the Higgs propagator and the decay channel. (Herein, we will adopt 
this factorisation in Narrow Width Approximation, thereby neglecting small corrections of 
0{Thi/Mhi), as Thi <^ for a light Higgs state.) In this process, effects of CP violation 
can occur through the aforementioned couplings in the production, through a possible mixing 
of Higgs states at one-loop and above in the propagator and through the same couplings in 
the decay. CP violation in the production of a Higgs state in the gluon-gluon fusion process 
at hadron colliders was studied first by [26] , choosing a parameter space region which is not 
sensitive to the CP mixing of the Higgs states, and later by [16, 28], including the presence 
of CP mixing of the Higgs states^. Effects of CP mixing in the propagator are discussed 
separately but in great detail in [20]. A detailed study of the other MSSM Higgs decay 
channels in presence of CP violation can be found in [15]-[20], [27]-[31]^. 

Results of a random parameter space scan to understand the general behaviour of the 
BR{Hi — > 77) for non-zero 0^ values are reported in [33]. It has been seen that about 
50% deviations are possible for M^^ around 104 GeV for 0^ = 100°, and an average of 30% 
deviation occurs over the mass range 90-130 GeV. Masses around and below 110 GeV show 
a decrease in the BR for non-zero 0^ values, while masses above this value show an increase. 
Certain individual parameter space points were also discussed in [33]. Specifically, \Af \ = 1.5 
TeV, \n\ = 1 TeV, tan/? = 20 was considered as a benchmark scenario. Then, by choosing 
(a) = 1 TeV and (b) M^^ = 250 GeV, it was demonstrated that the light stop ii has 
a strong impact on our decay mode, through the /i and the trilinear couplings Af, which is 
quite different from the effects due only to the (one-loop) change of the HiW^W~ coupling 
(discussed in [28]). In Ref. [33], it was also noticed that the effect of a light ti is in the 
opposite direction as compared to that due to modifications of the Hi coupling to the SM 
particles. 

In this article we will consolidate the results of [33] by a detailed discussion of the following 
points: 

^CP violation in vector boson associated production {VHi) is studied by [29]. 

^Where all these decay modes have been studied inclusively. A discussion of CP violation in exclusive 
four lepton final states via gauge boson decays can be found in Ref. [32] . 
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• The Higgs mixing matrix elements are discussed in detail showing the changes in mixing 
as the phase (p^ is varied. 

• The decay amplitudes due to individual (s)particles in the loops are presented in detail 
and the real and imaginary parts of the scalar and pseudo-scalar contributions are given 
separately. We will see that this unambiguously shows that a light stop contribution 
is comparable to the SM one, while the other sparticles have negligible effects even in 
scenarios when they are light. 

• Compared to [33], where only one particular value was considered for tan/3, |//| and 
\Af\, we present here results with four different tan/3 = 2, 5, 20, 50, two different = 
0.5, 1 TeV and two different \Af \ = 0.5, 1.5 TeV. 

We postpone the full analysis of gg/qq ^ Hi 'j'j to a forthcoming publication [34]. The 
outline of this paper is as follows. In Sect. 2 the CP mixing in the Higgs sector is explained 
in more detail. In Sect. 3 we analyse the phase dependence of Hi — > 77. We conclude in 
Sect. 4. 

2 Higgs Mixing in the CP-violating MSSM 

In the Higgs sector of the MSSM with explicit CP violation the CP- violating phases introduce 
non-vanishing off-diagonal mixing terms in the neutral Higgs mass matrix, which in the 
weak basis (0i,02,fl), where arc the CP-even states and a is the CP-odd state, may 
schematically be written as [6, 17, 27, 35] 

Here, Ai'g is a 2 x 2 matrix describing the transition between the CP-even states, Ai^ 
gives the mass of the CP-odd state while Mps = (A^|p)^ (a 1 x 2 matrix) describes the 
mixing between the CP-even and CP-odd states. The mixing matrix elements are typically 
proportional to 

MIp (X Im{nAf) (2) 

and dominated by loops involving the top squarks. As a result, the neutral Higgs bosons of 
the MSSM no longer carry any definite CP-parities. Rotation from the EW states to the 
mass eigenvalues, 

{(Pi,h,af^O{Hi,H2,Hsf, 

is now carried out by a 3 x 3 real orthogonal matrix O, such that 

O^MlO = diag(M^^, Mi^, M^.^) (3) 

with < < Mhs- As a consequence, it is now appropriate to parametrise the Higgs 
sector of the CP- violating MSSM in terms of the mass of the charged Higgs boson, Mh± , 
as the latter remains basically unaffected. (For a detailed formulation of the MSSM Higgs 
sector with exphcit CP violation, see Refs. [6, 27].) 
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Figure 1: Diagrams for Higgs decay into 77 pairs in the CP-violating MSSM: f = t,b,T; f = 

tl,2, &1,2, Ti,2- 

3 JTi ^ 77 in the CP-violating MSSM 

A Higgs boson in the MSSM decays at one-loop level into two photons through loops of 
fermions, sfermions, VF^ bosons, charged Higgs bosons and charginos, see Fig. 1. The 
analytical expressions for the respective amplitude along with relevant couplings in the CP- 
violating MSSM can be found in [36] and references therein. The amplitude has the form 

V Hi J 

where ki^2 are the momenta of the two photons and ei,2 are their polarisation vectors, which 
are conveniently written as e^j_ = — 2k!^{ks ■ er)/M'jj, (r 7^ s) and where (eie2/ci/c2) = 
^ij.vp(j^\^2^ik2- S] and are given by (retaining only the dominant loop contributions): 



SliMnd = 2 ^cQ}9fg'H,ff-:^F,finf) 

f=t,b,r,xt ,xi 



/j=tl,t2,6l,&2,Tl,T2 



-9HiWwFl{Tiw) — 9HiH+H--^Jj2~Fo{TiH±), 



V 



f=t,b,r,xt 'Xt 



For the expressions of the various couplings (gf's) and form factors {F{t)) we refer to [36]. 
Then the partial decay width is given by. 
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The decay mode Hi — * 77, i = 1,2,3, is discussed by Ref. [28] along with Higgs production 
through gluon-gluon fusion. However, that study was confined to MSSM parameter space 
regions with suitably heavy sparticles / and where CP-violating effects are only due 
to the changed SM particle (especially W^) couplings to the Hi and effects of sparticles 
in the triangle loops entering the decay amplitude are negligible. Here we examine the 
complementary region of MSSM parameter space with light sparticles so that they contribute 
substantially to the latter. In particular, we will show that, in the presence of non-trivial 
CP-violating phases, regions of MSSM parameter space exist where the couphngs of the 
Higgs bosons to all sparticles in the decay loops are strongly modified with respect to the CP- 
conserving MSSM, thereby inducing dramatic changes on the Hi 77 width and Branching 
Ratio (BR). 

We have analysed the Higgs decay widths and BRs with the publicly available Fortran 
code CPSuperH [36] version 2, which is based on the results obtained in Refs. [15]-[18] and 
the most recent renormalisation group improved effective-potential approach, which includes 
dominant higher-order logarithmic and threshold corrections, 6-quark Yukawa-coupling re- 
summation effects and Higgs boson pole-mass shifts [27, 37]. CPSuperH calculates the 
mass spectrum and decay widths of the neutral and charged Higgs bosons in the general CP- 
violating MSSM including the phases of ^4/ and jj,. Furthermore, it computes all the couplings 
of the neutral Higgs bosons i?i,2,3 and the charged Higgs boson H^ to SM particles and their 
superpartners. 

The open non-SM parameters of the model now include: the higgsino mass its phase 
(pn, the charged Higgs mass Mh±, the soft gaugino masses M^, the soft sfermion masses of the 
third generation M^q^ the (unified) soft trilinear coupling of the third generation 

1^4/1 and its phase (pAf 

In our analysis we fix the following SUSY parameters: 

• Ml = 100 GeV, M2 = 1 TeV, M3 = 1 TeV 

• = = = = MsusY = 1 TeV, 
whereas the following parameters are varied as given below: 

• tan/3 = 2, 5, 20, 50 

• \Af \ = 500 GeV, 1.5 TeV 

• 4>Af = 0° (whereas the CP-violating effects in the sparticlc sector depend on both 0^ 
and (pAf, the leading CP- violating effects on the Higgs sector, as stated above, are 
proportional to Xm{fiAf), and so we opted to fix (pAf to 0° and varied only 0^) 

• \n\ = 500 GeV, 1 TeV 

• 0^ = 0° - 180° 

• Mh+ = 100 - 300 GeV 
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• Mfj^ = 250 GeV (case with light ti) and M^j^ = 1 TeV (case with no hght sfermion) 

For this analysis, threshold corrections induced by the exchange of gluinos and charginos 
in the Higgs-quark-antiquark vertices [38, 39] were not included. While these corrections 
may change the actual values of the width and BR, we expect it to be the same for both 
CP-conserving and CP- violating cases. The situation will be different if 03, the phase of 
M3, could be non-zero. As mentioned in the introduction, we have considered the case of a 
common phase for the gaugino mass terms, which is rotated away making (a = 1, 2, 3) 
real. 

The mass of the lightest neutral Higgs boson, Mhi , is sensitive to the value of 0^ chosen. 
This dependence on 0^ along with that on other SUSY parameters is illustrated in Fig. 
2, where Mh^ is plotted against the charged Higgs mass for: tan/3 = 2, 5, 20, 50; = 
0.5, 1.5 TeV and = 1 TeV. Concerning the sparticles in the loop two cases are considered. 
The first case comprises a light m^^ ~ 200 GeV (corresponding to Mjy,^ = 250 GeV and 
Mq,^ = MsusY = 1 TeV) while all other sparticles are heavy, while in the other case m^^ is 
also taken in the TeV range with Mfj_^ = Mq^ = Msusy = 1 TeV. The top row in Fig. 2 
shows the sensitivity of Mh^ to tan/3 for |//| = 1 TeV and \Af \ — 1.5 TeV. While in the low 
tan/3 case the mass shift induced by the change in 0^ from 0° to 90° is about 10 %, in the 
case of tan /3 = 20 or above it is about 1 % or less. Notice that the relevant parameters here 
are \Af \ and and Mhj^ is found to increase with increasing \Af \ while its dependence on 

is basically negligible. In the bottom row of Fig. 2 we plot the Mh^ dependence on Mh± 
for two representative values of \Af\, 0.5 and 1.5 TeV, by keeping tan/3 = 20 and — 1 
TeV. 





0/. = O° 


0^ = 90° 


Mh+ (GeV) 


M„, (GeV) 


Mh, (GeV) 


M„, (GeV) 


Mh, (GeV) 


Mh, (GeV) 


Mhs (GeV) 


100 


40.7 


61.4 


120.7 


37.6 


61.2 


120.7 


120 


77.4 


90.3 


120.8 


75.8 


89.7 


121.0 


200 


120.5 


179.8 


183.8 


119.8 


176.7 


184.2 


250 


120.8 


232.1 


237.3 


120.2 


231.4 


237.6 


300 


120.9 


284.8 


289.6 


120.3 


284.2 


289.8 



Table 1: Selected values of Mh^ (i = 1, 2,3) for (p^ = 0° and (p^ = 90°. All SUSY parameters 
are as in Fig. 3 with = 250 GeV. 



The sudden shift in the dependence of Mh^ on Mh+ around = 150 GeV is under- 
stood in terms of the cross over in the mass eigenstates at that point. We have illustrated 
this in Fig. 3, where the masses of -f^i,2,3 are plotted against Mfj+ for tan/3 = 20, \Af \ = 1.5 
TeV and = 1 TeV, again with light and heavy stops. The cross over is a reflection of the 
changing compositions of the CP-indefinite mass eigenstates. Hi, H2, H3, with eigenvalues 
< < Mff^, in terms of the CP-definite gauge eigenstates, 0i, 02 and a. To explain 
this in a little more detail, let us denote the mass eigenvalues as mi, m2, m^, and the cor- 
responding eigenstates as hi, h2, /is, before ordering them from lightest to heaviest (i.e., hi 
need not be the lightest for all values of Mh+). Of the three mass eigenvalues, two (mi and 



7 



120 



100 



20, 



My^ = 250 GeV |n| = 1 TeV |A,| = 1 .5 TeV 
tanp = 20/ /_.,.,,^ tanf 



tan p = 2 



> 








a 


80 
















115 - 




60 




114 : 








113 J 




40 




93 i 








90 : 






/ /tan p = 50 


87 ': 



121.4 


- ' ■ ^ari ^ = so - 


121 




120.6 






100 



60 



40 



Mli^ = 1 TeV lul '=1 TeV |A,| = 1.BTeV 
tan p = 20 











tan p - 2 



' tan p = 50 



123.2 
122.4 


tan ^ = 50 - 










122.4 
121.6 




tan p = 20 - 



150 



200 
Mh* GeV 



250 



300 



■ 100 



250 



300 



90, 



My = 250 GeV |n| = 1 TeV tan p = 20 




1A,1 = 1.5 TeV 



|A,] = 500 GeV 



110 



100 



150 



J12Q. 



|A,| = 1.5TeV 



115 
110 
" 105 
100 
95 



250 



300 



90 



My^ = 1 TeV lul = 1 TeV tan p = 


20 






1 .5 TeV 






500 GeV 




^^^^^ ^ 


115.4 
115 


[A,l = 0.5 TeV 












114.6 










D o o o o o - 




122.4 








122 
121.6 


|A,|.1.5TeV 





100 



150 



200 
Mh» GeV 



250 



300 



Figure 2: Mass of the lightest neutral Higgs boson Hi against Mh+ for cf)^ = 0° (solid, red 
line) and (f)fj, = 90° (dashed, green line) with \Af\ = 1.5 TeV, \fi\ = 1 TeV and different 
values of tan P (top row) and tan/3 = 20, = 1 TeV and different values of \Af\ (bottom 
row), respectively. In the left column a light stop 200 GeV) is present for M^^ = 250 
GeV and Mq^ = Msusy = 1 TeV, whereas in the right column all sparticles are heavy (^ 1 
TeV) for = M^^ = Msusy = I TeV. 
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Figure 3: Mass of Hi 2^3 against Mh+ for tan/3 = 20 showing the cross over at Mh+ ~ 150 
GeV. Left column is for Mjy^ = 250 GeV (nii^ ~ 200 GeV), while the right one is for 
Mjj^ = 1 TeV (no light sparticle). Both plots are with \Af\ = 1.5 TeV, |/i| = 1 TeV. Red, 
blue and cyan curves represent (p^ = 0° , while green, magenta and orange curves correspond 
to 6,, = 90°. 
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Mh, (GeV) 


Mh, (GeV) 


Mh, (GeV) 


Mh, (GeV) 
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236.4 


237.2 
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289.2 
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122.2 
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Table 2: Same as Table 1, but with M^^ = 1 TeV. 



77^2) grow linearly with Mfj+ with almost the same slope, and lying close to each other. In 
the CP-conserving case one of these is a pseudo-scalar, while the other is a scalar. The other 
eigenvalue (ms) corresponds to a scalar eigenstate, and is more or less independent of Mh+. 
At Mh+ ~ 150 GeV all three CP-conserving eigenstates are degenerate with eigenvalues 
around 120 GeV. In the Mh+ ^ 150 GeV region hi is the lightest, while in the Mh+ ^ 150 
GeV region it is which is the lightest. When we order such that the lightest is Hi, there 
is a transition from Hi = hi to Hi = h^ around Mh+ = 150 GeV. For other values of tan /3, 
Af and the situation is very similar, with small shifts in the actual values of Mh+ and the 
degenerate mass where the cross over happens. In the CP- violating case with non-zero value 
of 0^, there is mixing between scalar and pseudo-scalar states. For larger values of M11+ the 
lightest Higgs state. Hi, is almost a pure scalar, hence we will not be subject to any CP 
violation through mixing. The only possible way to have CP violation here is through the 
Hif f* coupling, especially that of the stop quark. We restrict to regions of parameter space 
with small Mh+ where the effect of mixing as well as that due to a complex 0i,2//* coupling 
are present. At large Mh+ values there will be scalar /pseudo-scalar mixing in the heavier 
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Higgs states, H2 and H^. (We will not discuss the two heavier states in the present article 
though.) Tabs. 1 and 2 illustrate selected values of the Hi,H2 and if 3 masses for sample 
choices of Mh± when 0^ = 0° and 90° in the presence of a hght stop and otherwise, respec- 
tively. Mass shifts can typically be of a few percent, particularly for light Higgs masses, in 
the former case while they are negligible in the latter. 

Next we analyse the Higgs mixing matrix for the parameters tan/9 = 20, \Af \ = 1.5 TeV, 
lul = 1 TeV, as an example of a generic pattern over the entire MSSM parameter space. 
We show the mixing matrix elements in Fig. 4. In the CP-conserving case (0^ = 0°) Hi is 
mostly 4>i below Mh+ ~ 150 GeV and mostly 02 above it. H2 is the pseudo-scalar (a) below 
Mh+ ~ 150 GeV while above Mh+ ~ 150 GeV it is mostly 0i. H3 on the other hand is 
mostly 02 below Mh+ ~ 150 GeV and is the pseudo-scalar above this value. Indeed there is 
some mixing between 0i and 02 in the transition region. For the maximum value of 0^ = 90° 
the lightest (Hi) is mostly the pseudo-scalar below Mh+ ~ 150 GeV, H2 is mostly 0i and H3 
is mostly 02- Above this region Hi is mostly 02, H2 is mostly a and H^ is mostly 0i. There is 
of course some mixing (albeit small) between all the three states (especially in the transition 
region). For values of 0^ in between 0° and 90° mixing could be large, as demonstrated by 
the case of 0^ = 40° in Fig. 4. Similar features also happen between 0^ = 90° and 180° (the 
other CP-conserving value). Concentrating on Hi, the lightest eigenstate, we have plotted 
the relevant mixing matrix elements in Fig. 5 for two cases with and without the presence 
of a light ii, which shows that indeed the mixing is affected by the presence of a hght stop. 

These CP-mixing effects feed into the decay amplitude of Eq. (4) through couplings of 
the Hi^s to the SM and SUSY particles in the loop (see Fig. 1) at one-loop and tree-level, 
respectively. In Figs. 6-10 we show different contributions to the amplitude of Hi 77. 
Clearly, the SM contribution is dominant in all cases. Among the major contributions within 
the SM, that from the loop is about 5 times larger than the top quark contribution for 
the whole range of Mh+ (for the chosen set of SUSY parameters), while the bottom quark 
and tau lepton contributions are about an order of magnitude smaller over the lower range 
of Mh+ (except around 100 GeV) and negligibly small for larger values. Magnitudes of both 
the and top quark contributions grow with Mh+ , which is a reflection of the fact that 
these couplings are proportional to the mixing matrix element O21. When all SUSY states 
are heavy (Fig. 6), all sparticle (and if+) contributions to the real part of SJ are two or three 
orders of magnitude smaller than the SM term. The exception is the chargino contribution, 
which is at the most 10%. The effect of a light ti (Fig. 7) enters in two different ways. 
Firstly, it affects the SM couplings to Hi through loop corrections (compare top left plots 
in Figs. 6 and 7). This is more prominent when 0^ = 180°. Effects due to change in 0^ 
are different from the case with no light sparticle. Secondly, the ii contribution (top middle 
plot in Fig. 7) is now comparable (about 40% in the CP-conserving case) to that of the SM. 
Contributions of other sparticles are not changed much going from larger to smaller masses 
of the respective sparticle as shown in Fig. 8. In Fig. 9 contributions to the real part of Pi 
are plotted against Mh+. These terms come through the pseudo-scalar component in Hi. 
Since the squarks, sleptons and the charged Higgs boson do not couple to a, only SM objects 
and charginos contribute to . For the CP-conserving case it vanishes as expected. When 
0^ is non-zero Hi has an a component and there is a non- vanishing P^ , as illustrated by 
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the curve corresponding to 0^ = 90° in Fig. 9. Figure 10 shows the imaginary parts of Sj 
and which are sensitive to the value of 4>^. For the Hi mass range considered here only 
the SM contribution is complex. Again, being the contribution from a coupling to the 
(s)particles, its imaginary part vanishes in the CP-conserving cases. The Hi —>■ 77 width is 
not sensitive to the sign of Im P^, therefore the difference between the cases of light ti and 
no light sparticle is not very dramatic. 

From our MSSM parameter space scans, the following generic features on the sensitivity 
of width and BR to the CP-violating phase, 0^, have emerged. The effect is most pronounced 
around the crossing region {Mh+ ~ 150 GeV). This is expected since the scalar/pseudo-scalar 
mixing in Hi is largest here. For much higher values of Mh+ , Hi is purely a scalar. There 
can still be a CP-violating effect through its sfermion couplings though. The latter will be 
more visible when there is a light-sparticle in the loop (a stop, in particular). For Mh+ < 150 
GeV there is still sufficient mixing to have substantial difference in the BR. But in this region 
Mh^ is also changed by about 10 to 15% between 0^ = 0° and 0^ = 90°. We will however 
concentrate on the region Mh+ > 150 GeV, where Mh^ > 115 GeV for the parameter sets 
considered. Moreover, the effect of 0^ 7^ on Mh^ in this region is within 1 GeV, which is 
less than the experimental uncertainty expected at the LHC. In particular, we have learnt 
that the width and BR of the decay Hi — > 77 are very sensitive to the ti mass. The effect 
comes through both a modification of the HiW'^W~ couplings and the presence of a light 
ti in the triangle loops of the decay amplitude. 

As it is unfeasible to present all the results of our scan, we have picked out a few different 

discrete choices of the relevant MSSM parameters and plotted width and BR against Mh+ 
for the latter in Figs. 11-18. We choose five representative values between 0° and 180° for 
the phase of /i^ . Fig. 11 is with 1^4/1 = 1.5 TeV, = 1 TeV and tan/3 = 20. Comparing 
the two cases of M^^^ = 250 GeV and M^^^ = 1 TeV there is a qualitative difference in the 
sensitivity to 0^. This is also true for Fig. 12, which assumes \Af \ — 1.5 TeV and tan /3 = 20, 
but a smaller |/i| = 500 GeV. In the first case there is an increase in the BR over the region 
Mh+ > 150 GeV (and decrease over the region Mh+ < 150 GeV) as 0^ is switched on. This 
relative change with 0^ is maximised for some value of around 40°, beyond which the 
change in BR decreases again to about 50% at (f)^ = 180°. In the second case there is no such 
a trend as there is a 50% increase in the BR for 0^ = 90° at Mh+ ~ 200 GeV and the effect 
grows larger for 0^ > 90°. Other general features arc the following. The dependence on 0^ 
decreases with lower values of |^| as seen from Fig. 12. The value of 0^ with maximum BR 
in the presence of a light ti decreases compared to the case when |^| = 1 TeV. In contrast, a 
smaller value of \Af \ — 500 GeV (Figs. 13 and 14) keeps the picture qualitatively the same 
for the two cases of light and heavy ii. e.g., in the region Mh+ > 150 GeV, the BR increases 
with increasing 0^. But, while in the first case (M^^ = 250 GeV) there is a 50% increase 
for 0^ = 90° at Mh+ = 200 GeV, in the second case (M^^ = 1 TeV) there is a reduction of 
less than 20%. Again, the deviations can be substantially larger for 0^ > 90°. Sensitivity 

^To aid the reader, through multiple x-axcs and labels, we have included the mass of Hi in each plot 
corresponding to Mh± for the given choice of other MSSM parameters. Besides, in each figure we have 
zoomed near and above the aforementioned cross over point at Mh± ~ 150 GeV. Finally, as the BR is the 
measurable quantity, we dwell on this while only plotting the width for reference. 



11 



0,9 
0.8 




0.9 

08 


0.7 




07 


0,6 




06 


0=0,5 




("0,5 


0,4 




0,4 


13 




0,3 


02 
01 




0,2 
01 


n 









150 200 
M„.GeV 



200 
M„-GeV 







♦J-* 


0.5 












o" « 






-0.5 

.1 








Figure 4: Mixing matrix elements Oij vs. Mh+ (such that Hi 
with tan/? = 20, \Af \ = 1.5 TeV, = 1 TeV, Mq,^ = M^^ = 
TeV, Mfj^ = 250 GeV and (p^ as indicated in the plots. 



= On 01 + 02i 02 + a) 
= = MsusY = 1 



12 



















fii=* 






140 



























1 

0,9 
0.8 
07 
06 
J 0,5 
0,4 
03 
02 
0.1 



150 200 
M„-GeV 




200 

M„-GeV 

















fr™ 





















200 
M„-GeV 



200 
M„-GeV 





%< — 
(=9" 

9 =180° 












Figure 5: Mixing matrix elements On vs. Mh+ (such that Hi = On (pi + O21 02 + O31 a) 
for tan/3 = 20. Top row corresponds to Mjj,^ = 250 GeV while bottom row corresponds to 



M, 



Us 



1 TeV. All other parameters are as in Fig. 4- 




SoO 150 200 250 300 "°'^^100 150 200 250 300 '^■°°^100 150 200 250 300 

M^,* GeV M^^+ GeV Mh+ GeV 



Figure 6: Different contributions to Re (5*7) against the input parameter Mh+ with SUSY 
parameters as in Fig. 4 ond Mjy^ = 1 TeV. Top row: (from left) SM, ti^2 and bi^2- Bottom 
row: (from left) fi^2, X12 (^''^^ . 
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Figure 7: Different contributions to Re [SD against the input parameter Mh+ with SUSY 
parameters as in Fig. 4 o,nd Mjj^ = 250 GeV. Top row: (from left) SM, ti^2 o-nd hi 2- Bottom 
row: (from left) fi^2, ^'^'^ ■ 




Figure 8: Contributions to Re {Si) when the respective sparticle is light. Left plot is with 
Mf,^ =310 GeV (or m-^^ ~ 300 GeV), middle one is with M^^ = 310 GeV (or Mf, ~ 300 
GeV) and the right one is with M2 = 100 GeV (or M^± ~ 100 GeV). All other sparticles in 
the loop are of the order of TeV. The SM contribution in each case remains more or less the 
same as in Fig. 6 (no light sparticle). 
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of BR(iJi '-f'-f) to 0^ is however reduced considerably for lower values of tan j3, while the 
qualitative features remain the same, as is illustrated in Figs. 15-18 for tan (3 = 5. 

4 Conclusion 

In this paper, we have demonstrated that the decay channel Hi — > 77 is particularly suitable 
to probe the possible presence of CP-violating effects in the MSSM. This mode is in fact not 
only very sensitive to variations of the coupling HiW~^W~ pertaining to the dominant SM 
loop - induced from mixing amongst Higgs states via one-loop effects (as shown in previous 
hterature) - but also to contributions of a hght ii in the triangle-loop defining the decay 
process - which is in fact a tree-level effect induced by a complex // parameter (while the 
trilinear coupling At is taken real). 

In particular, our detailed analyses indicate that studies of the di-photon channel of 
a light Higgs boson (with mass below 130 GeV or so) found at the LHC may eventually 
enable one to disentangle the CP-violating case from the CP-conserving one, so long that 
the relevant SUSY parameters entering Hi — > 77 are measured elsewhere. This is not 
phenomenologically unconceivable, as the Hi — > 77 detection mode requires a very high 
luminosity, unlike the discovery of those sparticles (and the measurement of their masses 
and couplings) that enter the Higgs process studied here. Furthermore, while explicit CP 
violation could affect the mass of the hghtest Higgs state of the underlying SUSY model, 
we have restricted ourselves to regions of parameter space where - for an identical choice 
of all SUSY inputs but (^^, the only relevant phase in the scenarios considered here - the 
difference between the Mfji values in the CP-violating case and CP-conserving one are below 
the experimental uncertainty on the determination of such a quantity, so that it would not be 
possible to confirm or disprove the existence of complex parameters in the SUSY Lagrangian 
by solely isolating a Hi resonance. 

A complete analysis will eventually require to fold the decay process studied here in 
Narrow Width Approximation with propagator effects and the appropriate production mode 
(gluon-gluon fusion and Higgs-strahlung in this case) , where similar CP- violating effects may 
enter. This will be done in a future publication [34]. Finally, as argued in the Introduction, 
for the scenario we considered here, with very heavy squarks and vanishing trilinear couplings 
for the first and second generations, we can evade the EDM constraints on the CP-violating 
phases. A detailed analysis of this is also deferred to future work. 
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Figure 15: Similar to Fig. 11, hut with \ Af \ = 1.5 TeV, |/i| = 1 TeV and tan/? = 5. 
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Figure 16: Similar to Fig. 11, hut with \Af \ = 1.5 TeV, \fi\ = 0.5 TeV and tan/? = 5. 
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Figure 17: Similar to Fig. 11, hut with \ Af \ = 0.5 TeV, |/i| = 1 TeV and tan/? = 5. 
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Figure 18: Similar to Fig. 11, hut with \Af \ = 0.5 TeV, |/i| = 0.5 TeV and tan/5 = 5. 
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